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Model sets play a fundamental role in structure analysis of quasicrystals. The diffrac- 
tion diagram of a quasicrystal admits as symmetry group a finite group G, and there 
is a G-cluster C (union of orbits of G) such that the quasicrystal can be regarded 
£NJ \ as a quasiperiodic packing of interpenetrating copies of C. We present an algorithm 

which leads from any G-cluster C directly to a multi-component model set Q such 
that the arithmetic neighbours of any point x 6 Q are distributed on the sites 
of the translated copy x + C of C. Our mathematical algorithm may be useful in 
quasicrystal physics. 
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1 Introduction 



Model sets (also called cut-and-project sets) have been introduced by Y. Meyer 
[28] in his study of harmonious sets and later, in the course of the struc- 
ture analysis of quasicrystals, rediscovered in a variety of different schemes 
[17,18,20]. Extensive investigations [10,14,15,18,20,22,23,27,28,30,31] on the 
properties of these remarkable sets have been carried out by Y. Meyer, P. 
Kramer, M. Duneau, A. Katz, V. Elser, M. Baake, R.V. Moody, A. Hof, 
M. Schlottmann, J.C. Lagarias et. al. An extension of the notion of model 
set called multi- component model set, very useful in quasicrystal physics, has 
been introduced by Baake and Moody [2]. Model sets are generalizations of 
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lattices, and multi-component model sets are generalizations of lattices with 
colourings. 



Quasicrystals are materials with perfect long-range order, but with no three- 
dimensional translational periodicity. The structure analysis of quasicrystals 
on an atomic scale is a highly non-trivial task, and we are still far from a 
satisfactory solution. The electron microscopic images suggest the existence of 
some basic structural units which often overlap (interpenetrate), and of some 
glue atoms. The diffraction spectra contains sharp bright spots, indicative 
of long range order, called Bragg reflections. The reflections with intensity 
above a certain threshold form a discrete set admitting as symmetry group 
a finite non-crystallographic group G. In the case of quasicrystals with no 
translational periodicity this group is the icosahedral group Y and in the case 
of quasicrystals periodic along one direction (two-dimensional quasicrystals) 
G is one of the dihedral groups D$ (octagonal quasicrystals), Diq (decagonal 
quasicrystals) and Dyi (dodecagonal quasicrystals). 



The high resolution microscopic images of a quasicrystal with the symmetry 
group G show that we can regard the quasicrystal as a quasiperiodic packing of 
copies of a well-defined G-invariant finite set C (basic structural unit), most of 
them only partially occupied. From a mathematical point of view, C is a finite 
union of orbits of G, and we call it a G- cluster. In the literature on quasicrystals 
the term 'cluster' has several meanings [33]. Depending on the context, it may 
denote a structure motif (purely geometric pattern), a structural building 
block (perhaps with some physical justification), a quasi-unit cell [34] or a 
complex coordination polyhedron (with some chemical stability). In our case, 
C is a structure motif, perhaps without any physical justification. 



The purpose of the present paper is to present a mathematical algorithm 
which leads from any G-cluster C directly to a multi-component model set 
Q representing a quasiperiodic packing of interpenetrating copies of C, most 
of them only partially occupied. It shows how to embed the physical space 
into a superspace and how to choose a lattice L C in order to get by 
projection the desired local structure. Our algorithm, based on the strip pro- 
jection method and group theory, is a generalization of the model proposed by 
Katz and Duneau [18] and independently by Elser [10] for certain icosahedral 
quasicrystals. Since the multi-component model sets have several properties 
desirable from the physical point of view (they are uniformly discrete, rela- 
tively dense, have a well-defined density and are pure point diffractive) , our 
algorithm may be useful in quasicrystal physics. 
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2 Model sets and multi-component model sets 



In this section we review some definitions and results concerning the notions 
of model set and multi-component model set. 

Let E be a vector subspace of the usual /c-dimensional Euclidean space E& = 
(R k , (,)), where (x,y) = Yh=i %iVi and \ \x\ \ = \J(x,x), for anyx = (xj)i<j< fc , y = 
(yi)i<i<k- The subset B r (a) = {x G E \ \ \x — a\ \ < r} of E, where a G E, 
r G (0, oo), is the open ball of center a and radius r. 

Definition 1. Let A be a subset of E. 

(1) The set A is relatively dense in E if there is r G (0, oo) such that the ball 
B r (x) contains at least one point of A, for any x G E. 

(2) The set A is uniformly discrete in E if there is r G (0, oo) such that the 
ball B r (x) contains at most one point of A, for any x G E. 

(3) The set A is a Delone set in E if A is both relatively dense and uniformly 
discrete in E. 

(4) The set A is a lattice in E if it is both an additive subgroup of E and a 
Delone set in E. 

Definition 2. A cut and project scheme is a collection of spaces and mappings 



Ei <— ^- Ei © E 2 — ^> E 2 
U 
L 



formed by two subspaces Ei, E 2 of E fe , the corresponding natural projections 
tti, 7r 2 , and a lattice L in E x © E 2 such that: 

(a) 7Ti restricted to L is one-to-one; 

(b) 7r 2 (L) is dense in E 2 . 

Definition 3. A subset A of E n is a regular model set if there exist 

- a cut and project scheme (1), 

- an isometry X : E n — > E\ which allows to identify E n with Ei, 

- a set W 7^ satisfying the conditions: 

(i) W C E 2 is compact; 

(ii) W = mt(W); 

(iii) The boundary dW of W has Lebesgue measure 
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such that 



A = { -Kix | x e L, ir 2 x eW }. (2) 

By using the ^-mapping 

^(L) ^ E 2 : x^x* = n 2 ((n 1 )\ L y 1 x (3) 

we can re-write the definition of A as 

A = { x | x E tti(L), x* £ W }. (4) 

Model sets have strong regularity properties. 

Theorem 1. [30,31] Any regular model set A is a Delone set and has a well- 
defined density, that is, there exists the limit 

* (A ,? a V (5) 

where #(AnB r ) is the number of points of A lying in B r = B r (0), and vo\(B r ) 
is the volume of B r . 



In structure analysis of quasicrystals, the experimental diffraction image is 
compared with the diffraction image of the mathematical model A, regarded 
as a set of scatterers. In order to compute the diffraction image of the model 
set A, it is represented as a Borel measure in the form of a weighted Dirac 
comb 

u = 2 v( x )5x ( 6 ) 



where <p : A — > C is a bounded function and S x is the Dirac measure located 
at x, that is, 5 x (f) = f(x) for continuous functions /. In this way, atoms of 
quasicrystal are modeled by their positions and scattering strengths. In the 
case (the only considered in the sequel) when there is a function g : E 2 — ► C 
supported and continuous on W such that (p(x) = q(x*), that is, in the case 
^ = SxeA q(x*)S x one can prove [14,15,3,4] the following results: 

1. The measure ui is translation bounded, that is, there exist constants Ck so 
that 

sup \q(x)\<C k <oo (7) 

x eAn(t+K) 
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for all compact K C R™. 
2. The autocorrelation coefficients 




x, y G A n B, 



n 



q{x) g(y) 
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x — y = z 



exist for all z G A = A — A = {x — y \ x, y G A}. 

3. The set {z G A | rj(z) ^ 0} is uniformly discrete. 

4. The autocorrelation measure 



7c = v(z)S z 



(9) 



exists. 



The diffraction spectrum of A (the idealized mathematical interpretation of the 
diffraction pattern of a physical experiment) is related [14,15] to the Fourier 
transform % of the autocorrelation measure 7^ which can be decomposed as 



by the Lebesgue decomposition theorem. Here %(B) is the total intensity scat- 
tered into the volume B, (7 W ) PP is a pure point measure, which corresponds 
to the Bragg part of the diffraction spectrum, (7 w ) ac is absolutely continuous 
and (%) sc singular continuous with respect to Lebesgue measure. We say that 
A is pure point diffractive if 7^ = (%) pp , that is , if (%) sc = (%) ac = 0. We 
have the following result. 

Theorem 2 [14,31] Regular model sets are pure point diffractive. 

In the case of certain model sets used in quasicrystal physics as a mathematical 
model, the agreement between theoretic and experimental diffraction image is 
rather good [10,18]. 

The notion of model set admits the following generalization [2]. 

Definition 4. A subset A of E n is an m- component model set (also called a 
multi- component model set) if there exist 

- a cut and project scheme (1), 

- an isometry X : E n — > E\ which allows to identify E ra with E\, 

- a lattice M in E\ © E 2 containing L as a sublattice, 

- m cosets Lj = 9j + L of L in M, 
where j G {1, 2, m}, 



'lui ('lui)pp "I" if)iS)sc "I" ('111)) 



ac. 



(10) 
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- m sets Wj satisfying (i)-(iii), where j G {1,2, ...,m}, 



such that 



rn 



A = y {7rix | x e Lj, n 2 x e Wj}. 

3=1 



(11) 



Theorem 3. [2] Any multi- component model set is a Delone set, has a well- 
defined density and is pure point diffractive. 

The multi-component model sets have the property of finite local complexity, 
that is, there are only finitely many translational classes of clusters of A with 
any given size. The orbit of A under translation gives rise, via completion 
in the standard Radin- Wolff type topology, to a compact space X\, and one 
obtains a dynamical system (XA,M n ). The connection existing between the 
spectrum of this dynamical system and the diffraction measure allows one to 
use of some powerful spectral theorems in the study of multi-component model 
sets [24,25]. 



3 Model sets with predefined local structure 

Let {g : E n — > E n | g G G} be a faithful orthogonal M-irreducible represen- 
tation of a finite group G, and let 



be the G-cluster symmetric with respect to the origin generated by a finite set 
S C E n . For each g e G, there exist the numbers sf , sf, sf G { — 1; 1} and 
a permutation of the set {1, 2, k} denoted also by g such that 



Let e» = (eij)i<j< n , and let £j = (Sij)i<j< k , where 5^ = 1 for i = j, and = 
for % ^ j. 



C = [j Gx U |J = {ei, e k , -e u -e fe } 



(12) 



xes xes 



= s 9(i) e s(i) for a11 i e {I- 2 ? •••> 



(13) 



Lemma 1. [5,6] T"/ie formula gej = s 9 g ^e g ^) defines the orthogonal represen 
tation 




(14) 
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of G in Efc. T/ie subspace 

E = { (< «, ei >)i<i< fe | u G E n } (15) 

0/ Efc is G-invariant and the vectors 

Wj = « _1 (e^)i<i<fe j G {1, 2, n} 



where k = y (en) 2 + (e2i) 2 + ... + (e^i) 2 , /orm an orthonormal basis o/E. 

Lemma 2. [5,6] T/ie subduced representation of G in E is equivalent with 
the representation of G in E„, and £ne isomorphism of representations 

j :En ^ E Ju= (k- 1 ^,^))^.^ (16) 



mt/i the property I (a±, ot2, ot n ) = aiiui + ... + a n w„ allows us to identify the 
'physical' space E n win £ne subspace E o/Efc. 

6j Tne matrix of the orthogonal projector 7r:Efc — >E^ corresponding to E in 
ine frasz's {£1, £ 2 , £&} z' s 

7r =(«" 2 ^'^) 1 < i ,< fe 



Tae lattice L = /tZ fc C E fc is G-invariant, 7r(/c£j) = Ze i7 iaai is, 
7r(/c£j) = ej i/ we to£;e into consideration the identification X : E n — ► E, and 

tt(L) =Z ei + Ze 2 + ... + Ze fc . (18) 



Let V be a G-invariant subspace of E&. Since the representation of G in E& is 
orthogonal and (gx, y) = (gx, g{g~ 1 )y) = (x, g _1 y) the orthogonal complement 

V 1 - = {x G E fe I (x, y) = for all y G V} (19) 



of V is also a G-invariant subspace. The orthogonal projectors IT, Tl 1 - : E fc 
E fc corresponding to V and V 1 satisfy the relations 

non = n non^o no 5 = 5 on 
n^n^n 1 n x on = o n^^jon 1 



for any g G G. 
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Theorem 4. 

a) 7/n(L) is dense in V then L n V = {0}. 

&j 7/n(L) is discrete in V then L fl V contains a basis ofV. 

c) 7/n(L) is discrete in V then n ± (L) is a lattice in V -1 . 

Proof, a) Let us assume that there is z G L fl V, z ^ 0. For each y G L the 
solutions x = (xi,X2, ■ °f the equation 

Zi(&i - yi) + z 2 (x 2 - y 2 ) + ••• + z k (x k - y k ) = 

form the hyperplane i/^ orthogonal to z passing through y. The hyperplane if,, 
intersect the one-dimensional subspace Wz = {az \ a G M} at a point corre- 
sponding to a = (y, z)/\ \z\ | 2 . Since (y, G /t 2 Z, the minimal distance between 
two distinct hyperplanes of the family of parallel hyperplanes {H y \ y G L} con- 
taining n(L) is k 2 /\ \z\ |. The set n(L) which is contained in union H = U^gl Hy 
can not be dense in V. Each point of Wz — H belongs to V but can not be the 
limit of a sequence of points from n(L). 

b) In view of a well-known result [9] concerning lattices in subspaces of E fe , 
there exist Ai, X 2 , ■■■ , A s in L such that {LTAi, IIA2, IIA S } is a basis in V 
and II(L) = ZLTAi + ZIIA2 + ...ZI1A S . We extend {Ai, A 2 , A s } up to a basis 
{Ai, A2, Afc} of Efe by adding new vectors A s+ i, ... , A& from L. For each 
i G {s + 1, s + 2, k} there are an, ... , cti S G Z such that 

LTAi = ctjiilAi + aj2nA 2 + ... + «j s riA s 
that is, 

n(Aj - ajiAi - a i2 A2 - ••• - ct is \ s ) = 0. 

The linearly independent vectors 

Vi = \i- a a Ai - ... - a is \ s ie{s + l,s + 2, k} 

belonging to L form a basis in V L . Since the coordinates vn, Vi 2 ,...,Vik of each 
vector Vi belong to kL the space V which coincides to the space of all the 
solutions x = (xi,x 2 , ■■■,Xk) of the system of linear equations 

v a xi + v i2 x 2 + ... + v ik x k = i e {s + l,s + 2, ...k} 

contains s linearly independent vectors v±, v 2 , ... , v s from L. They form a 
basis of V. 
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c) The vectors v±, v 2 , ... , Vk from L form a basis of E^. Since 

for i G {l,2,...,s} 
Vi for i e {s + 1, s + 2, £;} 

and the change of basis matrix from {ke±, ke 2 , ■ ke^) to {v±, v 2 , ■ has 
rational entries it follows that the entries of II in the basis {kei, KEk] are 
rational. If q is the least common multiple of the denominators of the entries 
of IT- 1 then n-'-(L) is contained in the discrete set (n/q)Z k . ■ 

In order to obtain a description of the structure of Z- module n(L) we use the 
following result. 

Theorem 5. [9,32] Let : R k — > W be a surjective linear mapping, where 
I < k. Then there are subspaces Vi, V 2 ofM. 1 such that: 

a) R l = Vi © V 2 

b) </)(z k ) = <f)(z k ) n Vi + 0(z fe ) n v 2 , 

c) 0(Z fc ) fl V 2 is a lattice in V 2 , 

d) 0(Z fc ) fl Vi is a dense subgroup ojV\. 

The subspace Vi in this decomposition is uniquely determined. 

Theorem 6. IfV is a G-invariant subspace o/E^ and ifU is the corresponding 
orthogonal projector then there exist two subspaces Vi, V 2 such that: 

a) V = Vi © V 2 , 

b) n(L) = n(L) n Vi + n(L) n v 2 , 

c) II(L) fl V 2 is a lattice in V 2 , 

d) I1(L) fl Vi is a Z-module dense in V±. 

The subspace V\ is uniquely determined and G-invariant. 

Proof. The existence of the decomposition follows directly from the previous 
theorem. It remains only to prove the G-invariance of V\. For each x G Vi 
there is a sequence (£j)j>o in II(L) such that x = \im.j_ 00 C i j and ^- ^ x for 
all j. The transformation g : E& — > E^ corresponding to each g G G is an 
isometry and g(L) = L. Therefore 

gx = lim^oo g{,j g^ ^ gx 

ff (n(L)) = n^(L)) = n(L) 
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Fig. 1. Left: The decompositions E ( . = E©E 1 = E®E'©E" = £®E". Centre: A 
one-shell L>g-cluster C. Right: A fragment of the set Q defined by C. 

whence gx G V\. ■ 

In view of this result, there is a G-invariant subspace E' and a subspace V 
such that E- 1 = E' © V, 7r ± (L) n E' is a Z-module dense in E', and tt^L) n V 
is a lattice in V. Since E^ and E' are G-invariant, the orthogonal complement 
(see figure 1) 

E" = {i G E 1 (x, y) = for all y G E'} (20) 

of E' in E 1 - is also a G-invariant space. For each x G there exist 7r:r G E, 
G E' and x" G E" uniquely determined such that x = nx + x' + x" . The 
mappings 

7r' : E fc — >■ E fc : x i— >■ 7r'x = x' 

(21) 

tt" : E fc — ► E fc : a; i-> 7r"x = x" 

are the orthogonal projectors corresponding to E' and E". 
Theorem 7. The Z-module 7r(L) is either discrete or dense in E. 

Proof. In view of theorem 6, there is a G-invariant subspace Vi and a sub- 
space V2 such that E = Vi © V2, 7r(L) fl Vi is dense in Vi, and 7r(L) Pi V2 is a 
lattice in V2. Since the representation of G in E is irreducible we must have 
either Vi = {0} or Vi = E. ■ 

Let 7T = 7T + tt', £ = 7f(E fc ) = E © E', and let p : f — ► E, p' : 5 — >■ E' be 
the restrictions of n and n' to £. 

Theorem 8. The Z-module C = 7f(L) a lattice in £. 
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Proof. If ^(L) is discrete in E- 1 then E' = {0}, and in view of theorem 4 
the projection C = 7r(L) of L on £ = E is a lattice in £. 

If 7r J -(L) is not discrete in E 1 - then for any p G (0, oo) the dimension of the 
subspace V p generated by the set {k l x \ x G L, | |7r x a;j | < p} is greater than 
or equal to one. More than that, g < p' V p C V p ', and there is po G (0, oo) 
such that V p = V eo for any p < p . We have E' = V po . Since 



tt^o; £ E' 



|7r"x|| > po 



7r"(L) is a lattice in E". From theorem 4 it follows that £ is a lattice in 8. 
■ 

Theorem 9. Ifn(L) is dense in E then 
E £ E' 

U (22) 
is a cn£ and project scheme. 



Proof. From the relation p'(C) = tv'((tv + 7r')(L)) = 7r'(L) and the definition 
of E' it follows that p'(£) is dense in E'. 

Let x,y G £ with px = pi/, and let x, y G L be such that x = nx and y = ny. 
From p£ = p^ it follows nx = iry, whence x-i/GE 1 . Since L fl E" contains 
a basis of E" and L fl E' = {0} we must have x — y G E", whence x = y. 
Therefore p restricted to £ is one-to-one. ■ 



Let W = 7r J -(W) be the projection of the hypercube 

W = v + {(x t ) 1 < l < k \0<x t <K} (23) 

where the translation v G E' is such that no point of 7r J -(L) belongs to dW. 
Theorem 10. If 7r(L) is dense in E then the set 

Q = {nx | x eh, n L x G w} (24) 
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is a multi- component model set. 



Proof. Let 9 = {0< | % G Z} be a subset of L such that n"(V) = tt"(0) and 
7r"6i 7^ 7r"0j for i ^ j. The lattice L is contained in the union \J ie %£i of all the 
cosets Si = 9i + S of £ in The lattice L = L fl £ is a sublattice of £. Since 
L H £j = 0j + L the set Ci = 7r(L PI Si) = n:9i + L is a coset of L in C for any 
i £ Z. Since 7r"(L) is discrete in E", the intersection (see figure 1) 

w< = w n Si = w n tt ± (^) = tt^w n c tt"^ + e' 



is non-empty only for a finite number of cosets £j. By changing the indexa- 
tion of the elements of G if necessary, we can assume that the subset Wj = 
7r'(Wj) = 7r'(W H Si) of E' has a non-empty interior only for % e {1, ...,m}. 
The 'polyhedral' sets Wj satisfy the conditions (i)-(iii) from definition 3, and 

m 

Q = |J {nx \ xed, ir'x £ } . (25) 

i=l 



The set Q is a union of interpenetrating copies of the starting cluster C, most 
of them only partially occupied. For each point nx e Q the set of all the 
arithmetic neighbours of nx 



ny 



y G {x + KEi, ...,X+ KEk, X — KEi, X — KSk] 
7T-4/ G W 



is contained in the translated copy 

{nx + e 1 , nx + e^, nx — ei, nx — e^} = 7ra; + C 



of the G-cluster C. An algorithm and some software for generating such pat- 
terns is available via internet [8]. Some examples are presented in figures 1 and 
2. In each case, the quasiperiodic pattern is a packing of partially occupied 
copies of the corresponding cluster. One can remark that the occupation of 
these copies seem to be very low in the case of a multi-shell cluster. 

The number a is called a scaling factor of Q if there is y G E such that Q is 
invariant under the affine similarity [26] 

A : E — ► E Ax = y + a(x - y) (26) 
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Fig. 2. Left: A one-shell -Di2-cluster and a fragment of the corresponding quasiperi- 
odic set. Right: A fragment of the quasiperiodic set defined by a two-shell 
Dio-cluster. 

that is, if A(Q) C Q. In this case we say [26] that y is an inflation center 
corresponding to a, and A is a self- similarity of Q. The definition (25) offers 
some facilities [8,18] in the study of the self-similarities of Q. 

In view of the theorem 3, each quasiperiodic set defined by the above algorithm 
is a Delone set, has a well-defined density and is pure point diffractive. 



4 An example 

In order to illustrate the algorithm presented in the previuos section we 
consider the dihedral group D w = (a, b | a 10 = b 2 = (ab) 2 — e), the two- 
dimensional representation 

a(a, (3) = (a cos | — (3 sin |, asin | + (3 cos |) 
b(a,P) = (a,-P) 

and the -Di -cluster C = Dio(l, 0) generated by the set 5 = {(1, 0)}. 

The action of a and b on C generate the orthogonal representation of D w in 
E 5 

a(x ± , x 2 , x 3 , x 4 , x 5 ) = (-x 3 , -x 4 , -x 5 , -x 1 , -x 2 ) 
b(x 1 ,x 2 ,x 3 ,x 4 ,x 5 ) = (x 1 ,x 5 ,x 4 ,x 3 ,x 2 ). 
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The matrices of the projectors corresponding to E, E', E" are tt = M (|, — y, — 
n' = M (§, -£), n» = M (|, J, |), where r = (l + v^/2, t> = {l-y/E)/2 ' 



and 



/3 a /3 7 7 

7 P a (3 7 

7 7 /? a /3 

{0 1 1 P a J 



(28) 



In this case « = ^5/2, and the projection of L = /tZ 5 on the space 
£ = E © E' = {(xi, x 2 , x 5 ) | xi + x 2 + ... + a; 5 = 0} 



is the lattice £ = Z%i + Z%2 + Z% 3 + Z%4, where Xj = ( n + ^'X^')- ^ 
W = 7r- L (W) is the projection on E x = E' © E" of a hypercube 

W = v + {(xi, x 2 , x 3 , x±, x 5 ) | < Xj < k for any j} 



with dGE' chosen such that no point of ^(L) belongs to <9W then 

Q = {ttx I ieL, n L x e w} (29) 

is the set of all the vertices of a rhombic Penrose tiling [18]. 
The lattice L is contained in the union Ujez £? °f subspaces 

£?' = {(^1)^2, ...,^5) I ^1 + ^2 + ••• + %5 = jit} = 0j + £ 



where % = (jn, 0, 0, 0, 0) G JL. Since JL n £j = 9j + L, the set 
Cj = vr(L n £j) = n9j + L = jxi + L 



(30) 



is a coset of L = L (7 £ in C, for any j G Z. 

The set Wn£,- is non-empty only for j G {0, 1, 2, 3, 4, 5}, but W,- = 7r'(Wn£,-) 
has non-empty interior only for j G {1,2,3,4}. Let ? C E' be the set of all 
the points lying inside or on the boundary of the regular pentagon with the 
vertices tt'(k, 0, 0, 0, 0), 7r'(0, k, 0, 0, 0), 7r'(0, 0, 0, 0, k). One can remark that 
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Fig. 3. The cluster C = Dio(l, 0) and a fragment of the quasiperiodic pattern defined 
by this cluster. The positions of the points are indicated by using the symbols ., °, 
* and > in order to distinguish the four components of this model set. 

Wi = v + V, W 2 — v — tV, W 3 = v + tV, W 4 = v— V, and express the pattern 
Q as a multi-component model set (see figure 3) 

4 

Q= \J{px | are A, P'xe Wi}. (31) 

This definition is directly related to de Bruijn's definition [28]. 

The particular case of the icosahedral group, very important for quasicrystal 
physics, has been presented (with direct proofs) in [7]. 



5 Concluding remarks 

Quasicrystal structure analysis comprises the determination on an atomic scale 
of the short-range order (atomic arrangement inside the structural building 
unit) as well as the long-range order (the way the structural building units 
are arranged on the long scale). The rational approximants (periodic crystals 
with the same building unit as the considered quasicristal) provide a powerful 
way to determine the short-range order, but the description of the long-range 
order is still a major problem [33]. Only in the case of a few decagonal and 
icosahedral phases, the electron microscopic and diffraction data have allowed 
us to have an idea about the structure of the building unit and long-range 
ordering. 

The atomic structure can not be extracted directly from the experimental 
data. One has to postulate a structure and to compare the forecasts with the 
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electron microscopic and diffraction data. There exist several attempts in this 
direction: 

- Elser & Henley [11] and Audier & Guyot [1] have obtained models for icosa- 
hedral quasicrystals by decorating the Ammann rhombohedra occuring in 
a tiling of the 3D space defined by projection [10,18]. 

- In his quasi-unit cell picture Steinhardt [34] has shown (following an idea of 
Petra Gummelt [13]) that the atomic structure can be described entirely by 
using a single repeating cluster which overlaps (shares atoms with) neigh- 
bour clusters. The model is determined by the overlap rules and the atom 
decoration of the unit cell. 

- Some important models have been obtained by Yamamoto & Hiraga [35,36], 
Katz & Gratias [19], Gratias, Puyraimond and Quiquandon [12] by using 
the section method in a 6D superspace decorated with several polyhedra 
(acceptance domains). 

- Janot and de Boissieu [16] have shown that a model of icosahedral qua- 
sicrystal can be generated recursively by starting from a pseudo-Mackay 
cluster and using some inflation rules. 

Most information about the type of quasiperiodic long-range order is in the 
very weak reflections. The number of Bragg reflections we can observe is too 
small for an accurate structure description [33]. The experimental devices 
allow us to obtain diffraction patterns and to have a direct view on some 
small fragments of the quasicrystal. These data provide easy access to the 
symmetry group G and allows us to look for an adequate cluster C. The 
pattern Q obtained by using our algorithm is exactly defined and has the 
remarkable mathematical properties of the patterns obtained by projection. 
Each point of our pattern (without exception) is the center of a more or 
less occupied copy of C, but unfortunately, in the case of complex clusters 
the occupation is extremely low for most of the points of Q . Therefore, our 
discrete quasiperiodic sets can not be used directly in the description of atomic 
positions in quasicrystals. 

The cluster C can be regarded as a covering cluster, and Q as a quasiperiodic 
set which can be covered by partially occupied copies of a single cluster. This 
kind of covering is different from the covering of Penrose tiling by a deco- 
rated decagon [13] proposed by Gummelt in 1996 or the coverings of discrete 
quasiperiodic sets presented by Kramer, Gummelt, Gahler et. al. in [21]. In 
our case the generating cluster C is a finite set of points and the quasiperiodic 
pattern is obtained by projection. In [13,21] the covering clusters are congru- 
ent overlapping polytopes (with an asymetric decoration) and the structure 
is generated by imposing certain overlap rules which restrict the possible rel- 
ative positions and orientations of neighbouring clusters. When the theory 
from [13,21] is applied to quasicrystals, atomic positions are assigned to the 
covering clusters. 
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There are some indications that stable clusters are smaller than the basic struc- 
tural units seen on electron microscopic images, and it is believed that larger 
clusters automatically introduce disorder. The defects occuring in the tilings 
constructed from electron microscopic images show that a certain amount of 
disorder either in glue atoms or in the clusters seems to be unavoidable [12]. 
In the case of Gummelt's approach, the transition from perfect to random 
quasicrystalline order is obtained by passing to relaxed overlap rules [29] . The 
frequency of occurrence of fully occupied clusters in our quasiperiodic patterns 
can be increased by a certain relaxation in the use of strip projection method, 
but this leads to some defects. In order to correct these defects one has to 
elliminate some points from interpenetrating clusters if they become too close. 
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